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ABSTRACT 

If galactic halos contain stellar remnants, the infra-red flux from the remnant progen- 
itors would contribute to the opacity of multi-TeV 7-rays. The multi-TeV 7-ray horizon 
is established to be at a redshift z > 0.034 by the observation of the blazar Mkn501 . 
By requiring that the optical depth due to 77 — > e + e~ be less than one for a source at 
z = 0.034 we limit the cosmological density of stellar remnants, VL rm < (2 — 4) x lO -3 /^ 1 
(/170 is the Hubble constant in units of 70 km sec -1 Mpc -1 ) and thus strongly constrain 
stellar remnants as a cosmologically significant source of dark matter. 
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1. Introduction 

The constituents of the dark matter observed in 
galactic halos are, to date, unknown. Gravitational 
microlensing experiments, MACHO and EROS, have 
detected several unexplained events in the direc tion of 
the Magellani c Clou ds (Machos) (Alcock et al. 
Renault et al. 1997, Palanque-Dclabrouille -eta^ 
Their combined results limit Macho masses 
range (1(T 7 - 0.02) M^ to be less than 20% 
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Halo (Alcock et al. 1998). A standard interpreta- 
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tion of the microlensing results is that perhaps 30% 
of the galactic halo is composed of objects of mass 
roughly ~ 0.5M Q (Alcock et al. ), but this interpre- 
tation has several probl ems: Hubble Deep Fi eld sta r 
counts (Bahcall et al. |l996| ,Graff & Freese |l996a| ) 
an d an extrapolation of parallax data (Dahn et 
al. |1995| ) suggest that faint stars (O.O8M Q -O.2M ) 
and brown dwarfs (< 0.08 M ) contribute negligibly 
to the G alactic halo (<1% of the Halo mass) (Graff 



& Freese |l996b| ; Mera, Chabrier & Schaeffer |1996| ; 
Goul d, Flynn & Bahcall 1998). Gyuk, Evans & Gates 
998 have also argued that the Macho lenses cannot 
e a Halo or Spheroid population of brown dwarfs. 
Stellar remnants might make possible Macho can- 
didates; they have the right mass and are dark, but 
they also h ave their problems: Fields, Freese, and 



Graff (1998) extrapolated the Galactic population of 



Machos to cosmic scales and found a mass density 
that is comparable to the cosmic baryon density im- 
plied by Big Bang Nucleosynthesis; this result is in- 
dependent of the nature of the Machos. If the Ma- 
chos are stellar remnants, then the additional mass 
density of the gas left over from the progenitors be- 
comes a problem: virtually all of the baryons of the 
universe have to be recycled through the Machos and 
their progenitors. Over-pollution is difficult to avoid 
if these stellar remnants make up a significant fraction 
of the dark matter in galactic halos. It is problem- 
atic to hide the gas and/or metals ejected duri ng the 
formation of t he re mnants (Gibson and Mould 1997] : 
Fields et al. 1998 ). If carbon and nitr ogen d o not 
leave the stars, as suggested by Chabrier ( |1999| ), then 
the pollution is less problematic, although hel ium i s 
still very restrictive (Fields, Freese, and Graff |l999 ). 
These problems of baryonic mass budget and chemical 
overproduction are particularly severe for higher mass 
progenitors that give rise t o neutron stars (Venkate- 
san, Olinto & Truran |1999| ). 

In this Letter we constrain stellar remnant bary- 
onic dark matter by examining the infra-red radia- 
tion that is emitted during the evolution of the rem- 
nants' progenitors while on the main-sequence and 
during the subsequent red giant phase. This dif- 
fuse infra-red background (DIRB) contributes to the 



opacity of multi-TeV 7-rays allowing observations of 
TeV gamm a, ray sources to limit the DIRB ( Gould 
& Schreder 1966 , Ste cker. De Jager & Salamon 1992 , 
MacMinn & Primack |l996|) . HEGRA CT1 detections 



of mufti- TeV 7-ra ys from the blazar Mkn501 (Brad- 
bury et al. 1997) suggests that the universe is op- 
tically thin to 10 TeV 7-rays out to z — 0.034 and 
thus limits the DIRB (Funk et al. |1998| ). We use this 
multi-TeV 7-ray horizon to constrain the DIRB ex- 
pected from remnant halos. We show that the contri- 
bution of remnants to baryonic halo dark matter must 
be fi rm /i7 < 4 x 10 -3 where hro is the Hubble con- 
stant in units of 70 km s _1 Mpc _1 and f2 rm is the den- 
sity of stellar remnants in units of the critical density 
p c = 3H 2 /8ttG. Note that the existing constraints on 
the DIRB obtained by measuring the infrared back- 
ground directly with DIRBE (Hauser 1995), IRAS 
(Boulanger and Perault 1988), and other experiments 
are weaker than the constraints obtained by HEGRA 
for the energy range of interest. 

2. Measuring the DIRB with multi-Tev 7- 

rays 

The Universe appears to be optically thin to multi- 
TeV 7-rays out to a redshift of (at least) z ~ 0.03 
based on the detection of two blazars, Mkn421 (z = 
0.031) and Mkn501 (z = 0.034). Neglecting self- 
absorption, the dominant source of opacity for multi- 
TeV 7-rays (energy E(z) = (1 + z)E) is pair-creation 
off diffuse background photons (energy e(z) = (1 + 
z)e) which has a threshold of 
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where 6 is the relative photon scattering angle in the 
rest frame of the microwave background, -EtcV is the 
observed source photon energy measured in TeV, and 
we have fixed h = k = c = 1. The 77 — > e + e~ cross 
section is Bethe-Heitler and is plotted in Figure 1: 
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where 1 = 1 — (et/i/c) 2 and <tt = 6.65 x 10 25 cm 2 
is the Thomson cross section. 

As discussed below, we will be able to place lim- 
its on the DIRB in the energy range where the cross 
section is large, 0.03 — 3 eV. The optical depth r 77 at 
observed source energy E out to redshift z due to a 
comoving background spectrum (1 + z) 3 n(e)de is 
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Fig. 1. — The cross section for 77 - 
left to right, the three curves correspond to three dif- 
ferent Blazar photon energies, 10, 3 & 1 TeV. We can 
place limits on the DIRB where the cross section is 
large, in the range 0.03 — 3 eV. The blips around 10 
ev are numerical artifacts. 
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where fi = cos 9 and for a flat Universe, (l+z)dx/dz — 
H- 1 = H^lnfi + zf + n A }- 1/2 with H the Hubble 
parameter, (I the present matter density and Oa the 
present cosmological constant energy density in units 
of p c . Note that for sources at z=0.03, the value of 
f?A is irrelevant to the results. 

3. The DIRB Produced by Remnant Halos 

In order to determine the DIRB produced by the 
progenitor stars of remnant Machos, we start with 
stellar models of intermediate mass, (2 — 9)M Q , and 
low metallicity, Z = 10 -4 , generated specifically for 
this project. The models use the nuclear reaction 
cross-sections of Bahcall, Bau & Pinsonneault (1998). 
The equation of state is fully ionized in the interior, 
and the model uses the Saha equation at low tem- 
peratures (Guenther, Demarque, Kim & Pinsonneault 
1992). We used the model opacities of Alexander & 
Ferguson (1994) for T < 10, 000X, and OPAL opaci- 
ties for higher temperatures (Iglesias & Rogers 1996). 
We assumed grey atmospheres, as is appropriate for 
these high temperature models. Our models ran from 
zero age main sequence (ZAMS) to He core exhaus- 
tion. 

For comparison, Girardi et al. (1996) have also 
produced models for low Z isochrones, and find re- 
sults that are very similar; the similarity suggests 



that our models are theoretically robust. Girardi et 
al. claim that the post He core burning lifetime of the 
star is of order 0.3% of the main sequence lifetime, 
and can thus be ignored here. Our models contain no 
convective overshoot, and so conservatively underes- 
timate the total light emitted by stars. Were we to 
follow the convective overshoot system of Girardi et 
al. (1996), we would expect stars to live aproximately 
20% longer, emit a total of 20% more light, and thus 
our limits would be 20% more restrictive. 

In addition to the Z = 10~ 4 models discussed 
above, we also examined models with Z = 10 -8 
and with zero metallicity. We found no substantive 
changes in the stellar models, with the following ex- 
ception: for progenitors at the high mass end (9M©), 
stars of the extremely low metallicity of Z = 10~ 8 gave 
similar results to the low (Z = 10~ 4 ) models; how- 
ever, stars of strictly zero metallicity actually never 
evolved off the main sequence and produced about 
half as many infra-red photons as the low metallicity 
models. Although our limits would then be only half 
as severe, such a scenario of all of the stars having 
not even the merest whiff of metallicity is probably 
academic. 

Note that the mass range of our stellar models, 
(2 — 9)M Q , covers all the allowed white dwarf pro- 
genitors as well as some low mass neutron star pro- 
genitors. Remnant masses are take n from the models 
of Van den Hoek & Groenewegen ( 1997|) . Stars be- 
low 2M Q would leave remnants so bright that they 
would have been detected already (Graff, Laughlin & 
Freese 1998). The stellar models for stars above 9M Q 
are not considered as reliable as those for lower mass 
stars. In addition, as mentioned in the Introduction, 
Fields et al. (1998) showed that the mass budget and 
chemical abundance problems are particularly severe 
for the highest mass progenitors that give rise to neu- 
tron stars. 

We model the spectrum of light emitted by a star at 
a particular stage of stellar evolution as a black body 
spectral shape characterized by the effective temper- 
ature of the star. Unless there are features in the 
stellar photon spectrum that are much broader than 
the width of the cross section (as a function of infrared 
photon energy; the width ~ e t h), the optical depth is 
fairly insensitive to the spectral shape and is largely 
determined by the total number of photons with en- 
ergy above e t h (see Fig. 1). Any spectral features, 
even broad absorption bands, will have only a mini- 
mal effect on the optical depth. If anything, realistic 
spectra would show absorption at the ultraviolet end 
and re-emission of the same energy as more photons 
at the infrared end (see further discussion of this point 
below). Thus, a black body spectrum conservatively 
underestimates the number of photons produced by a 
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star []. 

Stars emit almost all their total energy during two 
distinct stages of evolution, main sequence (MS) and 
Helium core burning (HB). All other stages of stellar 
evolution are much too short lived or too dim to be 
significant in the total energy budget. Thus, we ap- 
proximate the total energy emitted by a star as the 
sum of two black bodies each marked by an average 
effective temperature (T) and a total emitted energy 
-Etot- These two quantities are determined by inte- 
grating over the stellar models as follows: 



energy e is 



Etot = I dtL{t) 

' stage 



(T) = E£ f dtL(t)T(t), 

J stage 



(4) 
(5) 



where T(t) is the effective temperature as a func- 
tion of time, and L(t) the luminosity. We summa- 
rize the stellar parameters adopted in Table 1, where 
we have given the initial mass of the progenitor, the 
remnant mass, the average effective temperatures of 
the main sequence (T)ms an d Helium core burning 
phases (T)hb, and the log of the corresponding total 
emitted energies of the two phases. 



Initial 


Remnant 


{T)ms 


log (.Ems) 


<r>HB 


log(£ H B) 


(M ) 


(M Q ) 


(°K) 


(ergs) 


(°K) 


(ergs) 


2 


0.68 


13260 


63.30 


5424 


62.65 


4 


0.91 


19840 


63.59 


12880 


63.22 


9 


1.84 


28940 


64.04 


19290 


63.49 



Table 1: Stellar parameters from our model 

The total number of photons emitted by a single 
star at energy e can thus be taken as the sum of the 
number produced by blackbody emission by the MS 
and RG phases, 



N(e) de 



Etots 

--{MS,RG} 
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{nh (T),) 4 exp(e/fc 6 (T) i ) - 1 



de. (6) 



Note that this equation has been normalized by re- 
quiring that / N(e)iede — E to t,i- We then integrate 
over the redshifts at which the radiation is emitted. 
The co-moving number density of white dwarfs is 
^rmPcninn- Thus we find that the co-moving num- 
ber density of background photons with present day 



1 An exception to this rule is low temperature zero metallicity 
stars which emit more energy at high frequencies than black 
body. However, our stars are at too high a temperature for 
this to be an issue. 



n(e) = dz Jl rm p c m rn JiV(e(l + z)) . 



(7) 



Results 



Since there is no observed break in the HEGRA 
gamma ray sp ectrum from Mkn 501 at redshift z=0.034 



(Funk et al. 1998), we can conservatively conclude 
that t 77 < 1 out to z=0.034. Using equations (p|) and 
(7j), we can then constrain the mass density ofwhite 
dwarfs due to the infrared emission by their progeni- 
tors. As mentioned previously, we consider a variety 
of progenitor masses in the range (2-9) Mq. Inspec- 
tion of Fig. 1 shows that we can obtain upper limits 
on the infrared photon density over the energy range 
3 x 1CP 2 — 3 eV, where the cross section for inter- 
action is sizable (this range would be extended for a 
hard spectrum). For each progenitor mass, we vary 
the gamma ray energy in the observed range (1-10) 
TeV to find the largest resultant value of optical depth 
r 77 . [Note that the gamma ray energy E is inversely 
proportional to the range of affected infrared photons 
(C.f. Eq. (1)).] 

We also consider a variety of redshifts at which 
the stellar light is emitted. Our limit in Eqn. (8) 
below applies to progenitor stars created at z < 60, 
but with some decreased sensitivity to stars in the 
mass range (3-4)M Q created at z > 15. Stars absorb 
photons above the Balmer cutoff at 0.5 eV, so that 
these photons are not available to scatter with the 
TeV gamma rays from Mkn 501. At z = 15, the 
Balmer cuttoff at 0.5 eV passes the threshold of a 
lOTeV photon at 0.03 eV; thus no light emerges from 
the star that can scatter with the TeV gamma rays. 
Since 2M Q stars emit most of their light just to the 
red of the Balmer cuttoff, we will not see the light 
from these stars emitted at z > 15. However, a low 
metallicity 2M© star will live for 600 Myr (assuming 
no convective overshoot), whereas the universe at z = 
15 is only 150/iyQ 1 Myr old. Thus, even if a 2M Q star is 
created at very high redshifts, it will emit most of its 
light after z = 15. A (4 — 9)M© mass star will be hot 
enough in its main sequence phase to emit much light 
at energies higher than the Balmer cutoff; we will see 
light from that star out to z=60, the redshift at which 
the Lyman cutoff is redshifted below 0.03 eV. There 
is still a small hole of decreased sensitivity to stars 
in the mass range 3 — 4M created at 15 < z < 60 
because of the absorption of light blueward of the 
Balmer cutoff; these stars are shorter lived than the 
2Mq stars and do emit a significant amount of their 
light within 150 Myr of their formation. We do not 
expect any significant amount of ba rvoni c objects to 
form before z = 60 (Tegmark et al. 1997 ). 
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Combining equations (|^) and (Q), we find a robust 
limit of 

n rm h 70 < (2 - 4) x 10~ 3 . (8) 

Since our limit depends on the number density 
of photons, which does not vary with redshift, and 
since we are sensitive to a large range of energies, 
our limit on il rm is not sensitive to the redshift at 
which the stellar light is emitted (except for the de- 
creased sensitivity discussed above for z > 15 and 
(3-4)M©). In addition, the limit is relatively insensi- 
tive to the star formation rate and to the initial mass 
of the progenitor star (with stellar mass between 2 
and 9 M©). The limits quoted above represent the 
full range of limits found for any star formation rate, 
initial mass function, and formation redshift within 
the broad bounds described above. The key variables 
that govern the severity of the limits are: higher mass 
stars give slightly more severe limits, and a "burst" 
star formation rate, which concentrates the emitted 
light within a narrow range of rcdshifts (8z/z < 1) 
will give more severe limits (for the burst, no photons 
redshift out of the testable energy range). 

One possible way to escape this limit is dust. It is 
possible that the progenitor stars are so enshrouded 
by dust that little of their light escapes, and if the dust 
is cool enough, that the light is re-radiated away at 
too long a wavelength to interact with the TeV 7-rays. 
Such a situation is most likely at high redshifts, where 
the dust could absorb the UV photons that later (by 
z=0.034) would have been redshifted into the infrared 
(IR). However, at lower redshifts (z < 10), dust could 
actually make the limits stronger: dust could absorb 
UV photons and reradiate them in the IR, causing 
much more absorption in the DIRB. Dust absorbtion 
is very model dependent, depending on the type of 
dust, the dust geometry and temperature, and beyond 
the scope of this work. 

5. Interpretation and comparison with other 
limits 

Our result constrains the cosmological abundance 
of stellar remnants. The microlensing experiments 
have only detected Machos around a single galaxy, 
the Milky Way. However, if we assume that the Milky 
Way is not a special galaxy, and that other, similar 
galaxies also have their coterie of Halo Machos, then 
the Universe should be filled with Machos. 



Dalcanton et al. (1994) searched directly for a cos- 
mological population of Machos by looking for a sig- 
nal of amplification of continuum emission of QSO's 
at high redshift. The fact that they did not find such 
an amplification allowed them to constrain all com- 
pact objects (not just remnants) in the mass range 

(O.1M - 1OM ) tO fi M acho < 0.1. 



Fields et al. (1998) examined the cosmic abun- 
dance of Machos. They found that a simple extrap- 
olation of the (supposed) Galactic population of Ma- 
chos to cosmic scales gives a cosmic density 



n 



Macho 



(0.0051 - 0.024) f gal h. 
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(9) 



Here f ga i is the fraction of galaxies that contain Ma- 
chos. To obtain an estimate of a lower limit to f ga i 
under reasonable assumptions, they considered the 
limiting case where only galaxies within one magni- 
tude of the Milky Way (i.e., Mmw ± lmag) contain 
Machos, and only spiral galaxies contain Machos, 
and found f ga i > 0.17. We consider these two results, 
the limit from Dalcanton et al. and Eq. (p[) , to be the 
most robust since they apply to all possible Machos, 
and not just stellar remnants. 

Fields, Freese & Graff (1998) also placed a much 
stronger limit on the cosmological density of stellar 
remnants by noting that the early universe had a low 
carbon and/or nitrogen enrichment. Using low metal- 
licity stellar yields of van den Hoek and Groenewegen 
(1997), they found that only about 10~ 2 of all baryons 
can have passed through the white dwarf progenitors, 
i.e., f2 rm < flprog < 10 _2 f2s < 4.5 x I0~ 4 h^ , where 
£Ib is the cosmic baryon density. However, as noted 
by Chabrier (1999), this limit depends on the carbon 
yields of zero metallicity intermediate mass stars. 

Our new limit in Eq. (^), based on the DIRB, 
is less general than limits on all Machos, but more 
robust than the carbon limit placed by Fields et 
al. (1998): although the carbon yields may be un- 
certain, intermediate mass stars certainly do produce 
light. Our new limit applies to remnants with any ini- 
tial mass function (between 2 and 9Mq) and any star 
formation rate (with a somewhat less restrictive limit 
for light emitted at 15 < z < 60 and for 2.5 - 4M ), 
and is thus extremely model independent. We can 
compare our new limit in Eq. (g) with the extrapola- 
tion of the Milky Way abundance of Machos placed 
by Fields et al. (1998) in Eq. (|). If all galaxies 
contain Machos in the same abundance as the Milky 
Way does, then the Machos cannot be stellar rem- 
nants (/170 = 1). Of course, it is impossible to place 
limits on the nature of the Milky Way Halo based only 
on cosmological limits: the Milky Way could be the 
only galaxy in the universe to have Machos. Still, 
based on this work, we can say that there must be 
some galaxies whose halos are not dominated by white 
dwarfs. 
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